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Abstract 

Autonomous multispecies systems witli more-tlian-two-neiglibor inter- 
actions are studied. Conditions necessary and sufficient for closedness of 
tlie evolution equations of the n-point functions are obtained. The average 
number of the particles at each site for one species and three-site inter- 
actions, and its generalization to the more-than-three-site interactions is 
explicitly obtained. Generalizations of the Glauber model in different di- 
rections, using generalized rates, generalized number of states at each site, 
and generalized number of interacting sites, are also investigated. 
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1 Introduction 



The principles of equilibrium statistical mechanics are well established. But, 
thermal equilibrium is a special case, and little is known about the proper- 
ties of systems not in equilibrium, for example about the relaxation toward the 
stationary state. Some interesting problems in non-equilibrium systems are non- 
equilibrium phase transitions described by phenomenological rate equations, and 
the way the system relaxes to its steady state. As mean-field techniques, gen- 
erally, do not give correct results for low-dimensional systems, people are moti- 
vated to study exactly-solvable stochastic models in low dimensions. Moreover, 
solving one-dimensional systems should in principle be easier. Exact results for 
some models on a one-dimensional lattice have been obtained, for example in 
[1-12]. Different methods have been used to study these models, including ana- 
lytical and asymptotic methods, mean field methods, and large-scale numerical 
methods. Systems with more than one species have also been studied [13-27]. 
Many of the arguments are based on simulation results. There are, however, 
some exact results as well. For most of the models studied, the interaction is 
between nearest neighbors. 

In a 10-parameter family of stochastic models with interactions between 
nearest neighbors has been studied. In these models, the A:-point equal time 
correlation functions {niUj ■ ■ ■ Uk) satisfy linear differential equations involving 
no higher-order correlations. We call these models autonomous^ in the sense that 
the evolution equations of ri-point functions are closed (contain only n- or less- 
point functions) . These linear equations for the average number of the particles 
{rii) have been solved. The same models were studied on lattices with boundaries 
in p9[ ] . It was shown that these models may exhibit dynamic and static phase 
transitions. The same idea has been generalized to multi-species models |22|] in 
one dimension with two-site interactions. There, conditions have been obtained 
that the Hamiltonian should satisfy in order that the evolution equation for 
correlation functions be closed. The set of equations for average densities can 
be written in terms of four matrices. These matrices are not determined uniquely 
from the Hamiltonian: there is a kind of gauge transformation one can apply 
on them which of course, does not change the evolution equation. A formal 
solution for the average densities of different species was found. The large-time 
behaviour of the average densities of different species was also studied. The time 
evolution equations for more-point functions, generally contain not only these 
four matrices, but also elements of the Hamiltonian, and to obtain a closed form 
for their solution is generally not easy. 

The Glauber dynamics was originally proposed to study the relaxation of the 
Ising model near equilibrium states. It was also shown that, there is a relation 
between the kinetic Ising model at zero temperature and the diffusion annihila- 
tion model in one dimension. There is an equivalence between domain walls in 
the Ising model and particles in the diffusion annihilation model. Kinetic gen- 
eralizations of the Ising model, for example the Glauber model or the Kawasaki 
model, are phenomenological models and have been studied extensively [30-41]. 

In this paper, autonomous multispecies systems with more-than-two-neighbor 
interactions are studied. Conditions necessary and sufhcient for closedness of the 
evolution equations of the n-point functions are obtained. As an example, we 
obtain explicitly the average number of the particles at each site for one species 
and three-site interactions. This is then generalized to the case where more than 
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throe site interact. As another c^xample, a generalization of the Glauber model 
is presented. In this generalization, the processes are the same as those of the 
ordinary Glauber model, but the rates depend on three free parameters, rather 
than one free parameter in the ordinary Glauber model. Finally this model is 
further generalized to the case where the number of interacting sites is more 
than three and the number of states at each site is more than two. 



2 Models leading to closed set of evolution equa- 
tions 

The models addressed are multi-species exclusion reaction-diffusion models. 

That is, each site is a vacancy or contains one particle. There are several 
kinds of particles, but at any time at most one kind can be present at each site. 
Throughout the paper, the system is assumed to be translationally invariant. 
Consider first a case whore the interaction is betwcx^n throe neighboring sites. 
Then the Hamiltonian describing the system can be written as 

L 

= X^-f^J,i+l,i+2- (1) 
i=l 

The number of sites is L and the number of possible states in a site is q (one of 
these states, for example the q-th one, may be the vacancy); different states of 
each site are denoted hy a, a = 1, • ■ ■ q. Introducing nf as the number operator 
of the particles of type a in the site i, we have 

1 



The average number of the particles of the typo a in the site i at the time t is 

K) = {SK\Pit)), (3) 

where \P{t)) := exp(tW) |P(0)) represents the state of the system at the time t, 

(5| = (s|0---0(s|, (4) 



and 



(11 •••!)• (5) 



So, the time evolution of (nf ) is given by 

^«) = (SKHTO). (6) 

The only terms of the Hamiltonian H entering the above equation are i/^.i+i i+i , 
Hi-i,i.i+i, and, Hi^u^i^i. The result of acting any matrix Q on the bra (s| 
is equal to that of acting the diagonal matrix Q on the same bra, provided 
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each diagonal element of the matrix Q is equal to the sum of all elements of 
the corresponding column in the matrix Q. So, the actions of (1 1 (g) n°')H, 
(1 (8) (g) l)H, and (n" (g) 1 (g) 1)H on (s| g) (s| g) (s| are equal to the actions of 
three diagonal matrices on (s| g) (s| (g) (s|. We use the symbol ^ to denote the 
equality of the action on (s| g) (s| (g (s|. We have 

(n" g) 1 g) 1)H - ^-^/37A n'' ® ® 
/37A 

(1 g) g) - J2 ^-^/37A g) (g 

/37A 

(lg)lg)n")i?-^ M;^^;,n''g)n^(gn^, (7) 
/37A 

where M^^;^'s are defined as 



1 Aa 
■^/37A 


q ffOlTUJ 


2 

-^/37A 




3 Aa 
-^/37A 


— q q TJ'T'^a 



An implicit summation (from 1 to q) over the a same subscript and superscript 
is always assumed. From these, (||) takes the form 

«) - '^^7A("f"7+i<2) + '^^7A«f-l"Z<i> + '^^7A("f-2«7-l'^'> (9) 

Generally, the right-hand side of (^) contains one-, two-, and three-point func- 
tions. (Note that n°"s are not independent.) We want to obtain a condition 
that only one point functions appear in the right-hand side. To do this, we 
consider the expression 

U = /a/37">f"fc> (10) 

and ask for the condition that the right-hand side is expressible in terms of 
linear combinations of n's, provided 

= 1. (11) 

It is obvious that if 



/q/37 — iFa + 2F 13 + sFf: (12) 

then the right hand side of ( |lo|) is expressible in terms of linear combinations of 
n's. To prove that this form for / is necessary as well, we just count the number 
of independent variables in f's satisfying the desired property. One can write 
nj in terms one 1 and other n"'s. Then it is seen that a general cubic form of 
71 's is expressible in terms of q'^ independent forms of n's, each containing no 
more than three n's. Of these, 1 -I- 3(g — 1) expressions (the monomials of degree 
zero and one of 3((?— 1) independent variables) are desirable. The coefficients of 
other monomials should be zero. So, from q'^ independent variables in /, there 
remains only 3{q — I) + 1 independent variables in /'s satisfying the desired 
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condition. It seems that the right-hand side of (|l^) contains more independent 
variables, 3q. But we note that the transformation 

iFa — > iFa+iB, (13) 

does not change the right-hand side of (|l^), provided 

J2^B^0. (14) 

i 

This means that there are 3 — 1 redundant variables in the expression (p^. So 
( p^ contains actually the correct number of independent variables, and hence 
is the most general form of / with the desired property. 

So, in order that (|^) be expressible in terms of only one-point functions, one 
must have 

It is noted that ]A's are not determined uniquely. Applying the gauge trans- 
formation 

]A'^ ^ = sp, with ^ (16) 

j 

does not change the right-hand side of (p^). 
If (^5|) is satisfied, (^ takes the form 

«)-(}^^+2^^+3^^)("f) 

+ lA${nl,)+lA'^p{nl,). (17) 

( p5|) in fact guarantees that the time-evolution equations of n-point functions 
contain only n- and less-point functions. In the simplest case, the one-species, 
each site is vacant or occupied by only one kind of particles. Then, the matrices 
^jA are two by two. 

One can do the same arguments for the case where more than three neigh- 
boring sites interact. Suppose the number of interacting sites is N. One defines 



'■^1-^.^= {U'^'^'j ■ (18) 

To ensure that in the time-evolution equation of one-point functions only one- 
point functions appear, one must have 

'^i-,.-T.)-^:- (19) 
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Here too, the gauge transformation ( jlq ) does not change the right-hand side 
of (p^), and hence the physics of the problem. It is to be noted that (|l9|) is 
also sufficient for n-point functions evolution equations to contain no more than 
n-point functions. 
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2.1 Some special cases 

We now consider some special cases. 

2.1.1 single species case 

In this case the matrices ^jA are two by two. The time-evolution equation for 
(rife) will then be 

(h,) = -a{n,) + f3{n,+i) + (3' {rn-i) + j{n^+2) + -,'{7^-2) + S. (20) 
5 can be eliminated using the redefinition 

Then, introducing the generating function 

00 

G{z,t) =Y,^i^\ (22) 

—00 

one arrives at, 

G{z, t) = (-a + fiz-^ + 0z + 70-2 + y^^) G{z, t), (23) 
the solution to which is 

G{z, t) = exp [t {-a + Pz'^ + P'z + jz'^ + j'z^)] G{z, 0). (24) 

Using 



^[z + z 



= E (25) 



fc= — 00 



exp 
one arrives at 

00 

x,{t) = e-"* Yl ^"-^-^^ {^t^) h (2t^) (/3'//3)('=-^-20/2(y/^)'/2^.(o). 
j,;=-oo 

(26) 

A similar procedure can be done for more-than-three-neighboring-sites inter- 
actions. The main difference will be the number of modified Bessel functions 
appearing in the expression. 

2.1.2 Generalizations of the Glauber model 

Consider a two-state three-neighbor interaction of the form 

U^,,= ^^^,, = 0. (27) 
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This means that, similar to the Glauber model, any site interacts only with its 
neighboring sites. The interactions are 



AAA - 


->A0A 


/"I 


000 - 


->0A0 




A%A - 


-^AAA 


Ai 


0A0 - 


-^000 


A2 


AA% - 


^A00 


Oil 


00A- 


^%AA 




AU - 


-^AA0 


Pi 


%AA - 


->00A 


02. 



(28) 

This is a generalization of the Glauber model. For the ordinary Glauber model 

J 



J 



Ml = = 1 — tanh 

Ai = At = 1 + tanh 
ai = "2 = A = /32 = 1 (29) 

The criterion ( |l9|) for the closedness of the evolution-equation for one-point 
functions results in the following relations between the rates in the generalized 
Glauber model. 

Hi - (3j - Xj for any i, j 

ai+f3i=a2+ fh (30) 

So there are four independent variables in terms of them the above eight pa- 
rameters can be expressed. One can write the expressions as 



Ml 


= A 


- B- 


C- 


c 


M2 


= A 


+ B 






Ai 


= A 


+ B + 




c 


A2 


= A 


- B 






ai 


^ A 


-B- 


c 




a2 


= A 


+ B + 


c 




/3i 


= A 


+ B + 


c 




f32 


= A 


-B- 


c. 





(31) 

Of course one of the parameters can be absorbed through a time-rescaling. The 
resulting evolution-equation for the average particle-number is 

(n,) = A + B - 2v4(?i,) + C(n,+i) + (32) 

which is easy to solve. 

This generalized Glauber model can further be generalized in two directions: 
when the number of the interacting sites is more than 3, and when the number 
of states of each site is more than 2. The first case means that the interaction 
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is in a block of length A'', resulting in the change of the state of a single specific 
site in that block. This rate of change depends on the states of this site and 
the states of the other — 1 sites. Let's label this specific site of the block by 
0. (Usually the length of the block is considered to be an odd integer {2k + 1), 
and the evolving site is assumed to be the central one.) Denote the state of the 
site i by (Ti, where (7^ can take the values 1 (particle) or (vacancy). Then the 
evolution equation for the average particle number is 

("o) ^ - ( H ^(1' ^) noll[l-n, + a,{2n, - 1)] \ 

\ i#0 / 

+ / ^ i?(0, a){l - no) n [1 - + ^^(2«. - 1)] ) • (33) 

Here the state of other interacting sites is denoted by a, R{cro, a) is the rate of 
change of the state of site 0, from ctq to 1 — (Tq: when the states of the other 
interacting sites is a. 

We are looking for those rates R{<Jo, (?) that make the right-hand side of this 
evolution equation a linear combination of (ni)'s. The claim is that the general 
form of these rates is 

R{ao, a)=A+ (-1)^« + ^ ^' '^^ ' ^^^^ 
Inserting this ansatz into the evolution equation ( |3^ ) we get 

{ho) = B + A{l-2no)+J2'^M)- (35) 

So it is clear that the ansatz (|^ leads into a closed set of evolution-equations 
for the average particle number. It remains to prove that this ansatz is the 
most general one satisfying this property. To see this, one considers (|3^). In 
the right-hand side of this equation there are 2^ terms (the expectation of 
monomials in terms of n^'s). Of these, we wish that the coefficients of all be 
zero, except for the coefficients of the constant term, and linear terms. So, there 
are 2^ — (A^-|- 1) equations to be satisfied for the rates (consisting themselves of 
2^ unknowns). This shows that the rates satisfying the desired property contain 
N +1 independent variables, and it is clear that the ansatz (Q) contains + 1 
independent variables. So it is the most general solution. 

Now consider the second generalization, when the number of possible states 
at each site is more than 2, say q, and each block consists of N sites. The state 
of the site i is denoted by Ci, which can take q values. That site the state of 
which is denoted by 0. The rate of change of the state of he site from /3 to a 
is denoted by Rp The evolution equation for the average numbers is then 



\ 3 t3=£a i^a I 

^(EE<.-on<'V (36) 
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Defining 



K,,:^-J2<B^ (37) 

f3^a 

the evolution equation reads 

«}=5]i?^,Vn^[]<A. (38) 

It is easy to see that for the right-hand side of this equation be expressible in 
terms of one-point functions, one should have 

^/3,ff = oRp + ^ iR^i (39) 
Note that this is nothing but 

"A'^p,, = j2W:,, (40) 



obtained in (|19[). In this case, other '^'s vanish. The evolution equation of the 
one-point functions is then 

(no")=5].iZ^(nf). (41) 

i 

Previous discussions shows that this is the most general form of the rates for 
them the evolution-equations of the one-point functions contain only one-point 
functions. Here too, there is a gauge freedom in choosing ii?'s, namely 

iR'^^ ,Rfj+,B''sp, with ^,5"-0. (42) 

i 

This is the same gauge freedom encountered with earlier, (p^). 
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